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ANOMALOUS  RESISTIVITY  ON  AURORAL  FIELD  LINES 


INTRODUCTION 

The  understanding  of  acceleration  processes  in  the  auroral  zones,  clue  to  parallel  electric 
fields  (  E  ),  is  closely  connected  with  the  problem  of  one-dimensional  anomalous  resistivity  and 
the  breakdown  of  runaway  acceleration.  The  one-dimensional  aspect  of  the  problem  comes  from 
the  fact  that  in  the  auroral  regions  of  interest  (i.e.,  h  >  1000  km),  D.e/cJe  >  1  (where  u>c,  £lc  are 
the  plasma  and  cyclotron  frequencies).  In  this  case,  the  Larmor  radius  of  the  electrons  is  less  than 
the  Debye  length,  so  that  even  in  the  case  of  ion  acoustic  turbulence  of  Debye  length  scales,  the 
magnetic  moment  of  the  electrons  (p  =  mufl2B)  is  an  adiabatic  invariant.  The  weak  magnetic  field 
case  (co,,^  IT,,)  results  related  to  anomalous  resistivity  are  not  directly  applicable.  This  is  obvious 
since  the  physical  process  by  which  anomalous  resistivity  results  when  coe  P-  is  the  interaction 
of  electrons  with  low-frequency  fluctuations  which  has  the  form  of  elastic  scattering.  The  rapid 
scattering  of  the  electrons  then  converts  the  directed  velocity  into  heating  and  the  appearance  of 
anomalous  resistivity.  For  £le  >  G3e,  the  elastic  scattering  is  not  allowed  since  p  is  constant  and  the 
electrons  can  only  slow  down  by  parallel  diffusion.  However,  as  noted  by  Petviashvilli  [1],  the 
parallel  diffusion  is  accompanied  by  the  formation  of  a  quasilinear  plateau,  which  reduces  the 
anomalous  friction  and  produces  trivial  resistivity  changes,  independently  of  the  amplitude  of  the 
low-frequency  waves.  A  corollary  puzzle  is  what,  if  anything,  inhibits  the  electrons  in  the  plateau 
region  as  well  as  the  negative  slope  region  from  freely  accelerating  [2,3] .  If  distributed  parallel 
electric  fields  are  to  exist  in  the  auroral  regions,  the  above  questions  should  be  answered. 

The  purpose  of  the  present  report  is  to  demonstrate  that  under  conditions  prevailing  in  the 
aurora,  the  conservation  of  the  first  adiabatic  invariance  is  violated  for  electrons  above  a  certain 
threshold  velocity.  This  effect  combined  with  trapping  a  major  portion  of  the  electron  distribution 
function,  by  large  amplitude  ion  cyclotron  waves  such  as  observed  by  S3-3  [4]  and  in  laboratory 
experiments  [5] ,  is  sufficient  to  answer  the  major  questions  posed  above  and  reproduce  most  of 
the  observed  features  of  the  auroral  energetic  electron  fluxes. 


PHYSICAL  MODEL 

Consider  a  plasma  in  a  uniform  magnetic  field  130  with  £2e  >  u>e,  and  an  electric  field  £  Q, 
parallel  to  Ba,  which  causes  the  electrons  to  acquire  a  drift  velocity  vd  large  enough  to  drive  an  ion 
cyclotron  wave.  There  are  two  basic  questions  to  be  answered.  First,  what  ultimately  limits  the 
growth?  Second,  what,  if  anything,  inhibits  the  electrons  from  freely  accelerating? 

Several  processes  have  been  discussed  in  the  literature  for  limiting  the  growth  [6,7] .  In  view  of 
the  recent  experimental  evidence  for  large-amplitude  (6 n/n  ^  1/2)  coherent  ion  cyclotron  waves  in 
the  aurora  [4]  we  consider  here  trapping  as  the  basic  stabilization  mechanism.  The  effect  of  large- 
amplitude,  low-frequency  density  fluctuations  on  the  plasma  resistivity  has  been  recently  examined 
by  Rowland  et  al.  [8] .  It  was  shown,  by  including  the  quasineutrality  effect,  that  most  of  the 
electrons  can  be  trapped  for  6n/n  >  1/4.  The  electron  distribution  function  is  composed  of  a  central 
part  carrying  little  or  no  current,  and  a  runaway  tail  of  density  nr  due  to  the  untrapped  electrons, 
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which  carry  the  current.  Therefore  J(t)  =  nrevr(t)  where  the  ratio  nrjn  is  controlled  by  the  level 
of  5n,/n;  as  shown  in  Rowland  et  al.  [8] ,  while  in  the  absence  of  a  mechanism  that  can  break  the 
adiabatic  invariance  of  the  energetic  electrons,  vr  *=  (eE0/m)t.  From  the  above  we  see  that  the 
existence  of  a  steady  or  quasisteady  state  reduces  to  the  answer  to  the  second  question,  i.e.,  the 
mechanism  that  can  inhibit  the  free  acceleration  of  runaways. 

To  answer  this  question,  we  examine  the  stability  of  distributions  which  are  composed  of  a 
cold  core  with  a  runaway  tail.  The  linear  stability  was  studied  by  Kadomtsev  and  Pogutse  (9] ,  who 
showed  that  low-frequency  oscillations  u>fe  =  u>ek  Jk  ^  where  is  the  wave  vector  along  the 
magnetic  field,  are  unstable  even  for  flat  distributions.  These  modes  are  principally  driven  by  the 
first  cyclotron  resonance  at  velocities  The  instability  has  a  growth  rate 

yk  *  nrln  tofc.  The  nonlinear  theory  of  the  instability  has  been  studied  extensively 

[3,10,11] .  It  was  shown  there  that  it  results  in  a  fast  isotropization  for  particles  with  velocities 
v  >vc  '  3( Qeltjue)vu,  (vu,  is  the  electron  thermal  velocity).  A  consequence  of  this  is  the  break¬ 
down  of  the  adiabatic  invariance  for  electrons  with  cu  >  i\. .  An  additional  slowing  down  process 
was  discussed  first  in  Papadopoulos  et  al.  [12]  and  confirmed  by  particle  simulations  by  Haber 
et  al.  [11] .  Since  only  particles  with  participate  in  the  resonant  scattering  and  have  their 

v  reduced,  electrons  will  tend  to  pile  up  at  vc  leading  to  the  formation  of  a  beam.  This  distribution 
function  can  suffer  friction  in  the  parallel  direction  due  to  wave  emission  not  only  by  the  cyclotron 
resonance  but  also  by  Cerenkov  type  beam  plasma  instability  at  the  lower  or  upper  hybrid  branch. 
Moreover  in  the  presence  of  a  constant  dc  electric  field  this  piling  up  can  be  enhanced  by  particles 
being  accelerated  up  from  lower  velocities.  This  reappearance  of  an  electron  beam  and  the  slowing 
down  of  the  electrons  in  the  parallel  direction  will  further  isotropize  the  electron  distribution 
function. 

Another  interesting  effect  that  is  seen  in  the  simulations  is  the  acceleration  of  particles  in  the 
opposite  direction  to  the  electric  field.  The  electrostatic  waves  that  are  generated  in  the  parallel 
direction  can  be  backscattered  by  the  ion  fluctuations  and  create  superthermal  tails  by  Landau 
damping.  These  particles  can  also  interact  via  the  normal  Doppler  resonance  [e  =  (u>fe  -  J2e)//cu] 
with  the  waves  that  are  pitch  angle  scattering  the  runaways  and  isotropize  the  counter  streaming 
particles. 


SIMULATIONS 

We  report  next  the  results  from  a  series  of  particle  simulations  that  combine  the  effects  of 
finite  ion  cavities  and  pitch  angle  scattering  at  the  anomalous  Doppler  resonance  (ADR).  The 
parameters  for  the  simulations  are  the  same  as  those  used  by  Haber  et  al.  [11  ] .  However,  instead  of 
starting  with  a  runaway  distribution,  we  start  with  a  Maxwellian  distribution  and  a  fixed  ion  cavity 
with  an  effective  depth  of  0.3.  We  will  apply  a  constant  dc  field  to  show  the  formation  of  the  run¬ 
away  distribution  such  as  was  seen  in  the  earlier  strictly  one-dimensional  simulations  of  Rowland 
et  al.  [8] .  We  will  then  continue  to  drive  the  electrons  to  determine  the  effect  of  the  pitch  angle 
scattering  on  the  further  acceleration  of  the  electrons. 

Figure  1(a)  shows  the  initial  electron  distribution.  The  magnetic  field  is  at  55°  to  the  .\  axis 
and  Qe  =  2oj(,.  The  electric  field  is  parallel  to  B  and  has  a  magnitude  of  0.01(m/e)e,(,()oJ(,  (where 
m(e)  is  the  electron  mass  (charge),  vleo  is  initial  electron  thermal  velocity,  and  u>(,  is  the  plasma 
frequency] .  Figure  1(b)  shows  the  electron  velocity  distribution  at  2000  oj~, 1 .  The  ion  cavity  has 
prevented  the  hulk  of  electrons  from  being  accelerated,  and  one  can  see  the  cold  dense  core  of  elec¬ 
trons  at  zero  velocity.  The  untrapped  electrons  accelerated  along  the  magnetic  field  forming  a  flat 
runaway  tail.  One  can  see  that  the  temperature  of  the  electrons  transverse  to  the  field  did  not 
increase,  and  up  to  this  point  the  acceleration  was  ID  parallel  to  the  magnetic  field.  However,  the 
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Kij*.  1  -  The  initial  \l)  acceleration  of  the  tart  runs  parallel  to  the  mat* 
netic  field  H  and  E„  is  at  55  to  the  axis,  (a)  initial  distribution,  (b)  distri 
bution  at  2000  u ;t,  1 .  The  ion  cavities  prevent  the  acceleration  of  the  bulk 
of  the  electrons. 


high-velocity  edge  of  the  distribution  has  become  greater  than  rc,  and  these  particles  started  to  pitch 
angle  scatter  as  seen  by  the  spreading  of  the  upper  edge  of  the  distribution. 

At  this  point  in  time,  that  2D  effects  started  taking  place,  the  simulation  was  stopped.  This 
state  was  used  as  an  initial  state  for  a  series  of  simulations  using  different  values  for  the  external  dc 
electric  field.  Our  earlier  theoretical  and  computational  work  |  Rowland  et  al.  |8]  I  demonstrated 
that  during  the  initial  ID  stage,  changing  the  strength  of  the  de  electric  field  did  not  affect  the 
formation  of  the  runaway  distribution  beyond  changes  in  the  timescale;  namely  with  a  stronger 
electric  field  the  runaway  distribution  formed  faster.  The  simulation  presented  in  detail  here  had  an 
electric  field  starting  from  2000  or’  1  of  magnitude  2  •  10“a(m(,/<})  i'((.()or(,.  As  the  pitch  angle 
scattering  continues,  the  formation  of  a  beam  in  the  parallel  direction  similar  to  the  one  seen  by 
Haber  et  al.  [11]  is  observed.  The  region  of  positive  slope  is  unstable  to  Cherenkov  interactions 
with  both  lower  and  upper  hybrid  waves.  This  acts  to  reduce  the  parallel  current  and  to  further 
symmetrize  the  distribution.  Figure  2(a)  shows  the  distribution  at  3800  or(.  1 .  No  further  parallel 
acceleration  of  the  electrons  is  observed.  Electrons  with  rn  >  i\.,  which  could  be  pitch  angle  scat¬ 
tered  to  large  i\  ,  were  slowed  down  by  the  beam  instability.  Thus  one  sees  electrons  with  r,  •  - 1\. 
but  with  t’j  >  This  isotropization  continues.  Figure  2(b|  shows  the  electron  distribution 
function  at  8800  u)~  It  is  composed  of  a  cold  dense  core  of  electrons  surrounded  by  a  hot, 
isotropic  cloud. 

Figure  3  shows  the  parallel  electron  distribution  at  2100.  3200.  1800.  8800  u>(.  1  (a.b.c.d). 

The  high-velocity  beam  seen  at  2400  u>(,  1  is  due  to  the  pitch  angle  scattering.  By  3200  u>(,  1  beam 
instabilities  have  flattened  the  runaway  distribution.  It  is  clearly  seen  that  the  instability  al  the 
anomalous  Doppler  resonance  prevents  electron  runaway.  If  the  high-velocity  electrons  had 
continued  to  freely  accelerate,  the  upper  edge  of  the  distribution  would  have  been  at  the  right- 
hand  edge  of  Fig.  3(d). 

Figure  1  shows  the  growth  of  the  parallel  drift  velocity  normalized  to  the  initial  electron 
thermal  velocity.  Note  that  the  time  axis  between  0  and  2000  oj(.  !  is  stretched  by  a  factor  of  five  to 
compensate  for  F.tl  being  five  times  larger  during  that  time.  Thus  the  time  axis  is  linear  in  terms  of 
the  effective  acceleration  time,  Knt.  Between  0  and  150  uj,.  1  the  ion  cavity  was  adiabatically  formed 
in  the  plasma.  At  150  u)(.  1  tin*  dc  electric  field  was  turned  on.  The  solid  line  shows  for  comparison 
the  rate  of  acceleration  of  the  electrons  in  the  absence  of  the  cavity.  Between  150  and  2000  uj(  1  we 
have  essentially  a  ID  system  parallel  to  the  magnetic  field.  The  cavities  prevent  the  acceleration  of 
the  bulk  of  the  plasma,  and  the  current  is  carried  by  a  small  fraction  of  high-velocity  electrons 
streaming  along  (lie  field  1  i nos.  At  approximately  2000  u)r  1  the  velocity  of  the  fastest  particles 
exceeds  and  pitch  ancle  scattering  begins;  the  high-velocity  electrons  are  heated  in  the  transverse 
direction  with  Cj  becoming  greater  than  rtl  and  the  adiabatic  invariance  is  broken.  The  current 
continues  to  increase  until  approximately  2300  ic(  '  when  the  instability  at  the  ADR  has  grown  to 
such  a  level  that  it  removes  energy  from  the  parallel  motion  and  slows  the  electrons  down  faster 
than  F(l  accelerates  them.  The  current  is  clamped,  and  strong  transverse  heating  is  observed.  Finally, 
at  much  longer  times  (  J  7000  u>(,  1 )  the  current  begins  to  increase  again  but  at  a  slower  rate  than 
during  the  ID  stage.  For  r  >  i\.  the  instability  heated  the  plasma  so  that  <)/7 <>t’ H  *  ()/7<h'| 
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Fig.  2  —  As  the  electrons  continue  to  accelerate  v  >  vc,  the  pitch  angle 
scattering  breaks  the  invariance  of  0  and  strong  transverse  heating  sets 
in.  (a)  t  =  3800  *,  (b)  t  =  8800  *. 


Fig.  3  —  The  parallel  electrons  distribution.  The  pitch  angle  scattering  prevents  the 
free  streaming  acceleration  of  the  electrons  t  =  2400,  3200,  4800,  8800  co" 1  (a,b,c,d). 


Fig.  4  —  The  average  electron  drift  velocity.  When  r  >  cr. 
the  pitch  angle  scattering  clamps  the  current. 
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(or  Tu  *  ).  The  electric  field  can  begin  to  increase  en  for  these  particles  but  the  rate  of  increase 

is  slowed  due  to  p'tch  angle  scattering  which  acts  to  keep  T  *  7’i.  The  increase  in  the  current  can 
be  seen  in  Fig.  3.  The  cutoff  at  vc  is  still  being  maintained  but  a  few  particles  are  accelerated  to 
higher  velocities. 

Figure  5  shows  the  log  of  the  electron  distribution  at  t  =  8800  u>~  1 .  A  cold  dense  core  of 
trapped  electrons  surrounded  by  the  hot  isotropic  accelerated  electrons  is  seen.  This  distribution 
has  marked  similarity  to  the  auroral  electron  distributions  measured  by  Kaufmann  et  al.  [13] . 


CONCLUSIONS 

The  acceleration  in  the  auroral  zones  is,  of  course,  a  spatial  problem.  While  temporal  simula¬ 
tions  such  as  reported  here  are  very  important  for  understanding  the  basic  microphysics  that  takes 
place  in  the  auroral  environment,  large  scale  macrophysics  simulations  are  needed.  The  present 
results  are  a  necessary  input  before  a  complete  understanding  of  the  total  system  dynamics  can  be 
gained.  Development  of  such  a  capability  is  presently  under  way.  We  present  a  qualitative  picture 
of  the  macrostructure  expected  on  the  basis  of  our  simulations.  The  first  point  is  that  for  the  pitch 
angle  scattering  to  be  effective  requires  £2e/ue  >  1.  Thus  one  would  expect  to  see  strong  dc  resis¬ 
tivity  and  parallel  electric  fields  between  2000  to  12  000  km.  The  critical  value  of  the  current  at 
which  strong  resistivity  appears,  increases  with  Z2e/u>e.  Thus  there  exists  a  maximum  critical  current 
that  can  be  carried  along  the  field  line.  Below  this  level  the  ion  cavities  prevent  the  bulk  of  the 
electrons  from  being  accelerated  and  tails  of  field  aligned  electrons  are  formed.  However,  when  this 
critical  current  is  exceeded  one  expects  the  emission  of  both  upper  and  lower  hybrid  waves,  large 
parallel  electric  fields,  and  a  hot  isotropized  electron  distribution  with  a  cold  dense  core.  The 
strongest  resistivity  and  hence  the  largest  electric  fields  will  appear  in  the  region  where  the  current 
is  clamped.  One  can  make  a  rough  estimate  of  the  initial  extent  of  this  region  by  assuming  Sr  =  vt 
where  r  is  the  speed  of  the  high-velocity  particles  and  r  is  the  time  over  which  the  current  is 
clamped.  For  the  simulation  shown  r  *  4.10;i  uj”  '.  Assuming  np  ==100  cm  and  a  velocity  of  a 
koV  electron  this  leads  to  a  distance  of  200  km.  This  is  of  course  only  a  rough  order  of  magnitude 
estimate  but  the  main  point  is  that  the  resistivity  and  the  parallel  electric  fields  should  extend  at 
least  initially  over  a  large  region  in  comparison  to  the  Debye  length  ( ==3  m).  This  is  in  agreement 
with  a  recent  analysis  of  satellite  data  of  potential  drops  along  auroral  field  lines  [14] .  The  dynamic 
spatial  behavior  is  a  subject  for  further  study  but  this  region  could  shrink.  For  the  simulation 
shown  in  detail  the  cu  electric  field  E0  =  2.7  •  10" 4  T}^nV2  (where  Tt,  is  in  eV,  np  is  plasma 
density  in  cm4,  and  Elt  is  in  V/m)  for  t  >  2000  to”  U  For  Te  =  20  eV  and  np  =  100,  Eit  —  10”  -  V/m 
and  ■>  E()  •  Sr  =  2  keV.  Based  upon  our  series  of  simulation  t  F.~  1 .  For  constant  S<p,  reduction 
of  Sr  would  lead  to  a  larger  t'() .  This,  in  turn,  would  lead  to  a  shorter  r  and  further  reduction  of  Sr. 
On  the  other  hand,  if  the  system  responds  by  raising  0  in  order  to  drive  a;  > ]cr.  an  increased  Elt 
and  a  shorter  Sr  will  result.  Such  questions  can  only  be  studied  with  a  large  scale  transport  code. 


Fig.  5  —  A  3/t  plot  of  (he  log  of  (hi’  electron  distribution  seen  in  Fig.  2(b). 
Note  the  cold  central  core  and  the  hot  isotropi/.ed  electrons. 
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Another  spatial  effect  presently  under  study  is  the  anomalous  transport  of  tin  hot  electrons  across 
the  magnetic  field  and  out  of  tin1  acceleration  region.  If  this  takes  place  at  a  fast  enough  rate  such 
that  7\  remains  less  than  T|j,  the  current  will  remain  damped.  In  this  case  t  will  be  determined  by 
the  timescale  for  this  perpendicular  transport. 
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